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Abstract. We address the basic postulates of quantum mechanics and 
point out that they are formulated for a closed isolated system. Since 
we are mostly dealing with systems that interact or have interacted 
with the rest of the universe one may wonder whether a suitable modi- 
fication is needed, or in order. This is indeed the case and this tutorial 
is devoted to review the modern tools of quantum mechanics, which are 
suitable to describe states, measurements, and operations of realistic, 
not isolated, systems. We underline the central role of the Born rule and 
and illustrate how the notion of density operator naturally emerges, to- 
gether the concept of purification of a mixed state. In reexamining the 
postulates of standard quantum measurement theory, we investigate 
how they may formally generalized, going beyond the description in 
terms of selfadjoint operators and projective measurements, and how 
this leads to the introduction of generalized measurements, probability 
operator-valued measures (POVMs) and detection operators. We then 
state and prove the Naimark theorem, which elucidates the connec- 
tions between generalized and standard measurements and illustrates 
how a generalized measurement may be physically implemented. The 
"impossibility" of a joint measurement of two non commuting observ- 
ables is revisited and its canonical implementation as a generalized 
measurement is described in some details. The notion of generalized 
measurement is also used to point out the heuristic nature of the so- 
called Heisenberg principle. Finally, we address the basic properties, 
usually captured by the request of unitarity, that a map transforming 
quantum states into quantum states should satisfy to be physically ad- 
missible, and introduce the notion of complete positivity (CP). We then 
state and prove the Stinespring/Kraus-Choi-Sudarshan dilation theo- 
rem and elucidate the connections between the CP-maps description 
of quantum operations, together with their operator-sum representa- 
tion, and the customary unitary description of quantum evolution. We 
also address transposition as an example of positive map which is not 
completely positive, and provide some examples of generalized mea- 
surements and quantum operations. 
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1 Introduction 

Quantum information science is a novel discipline which addresses how quantum 
systems may be exploited to improve the processing, transmission, and storage of 
information. This field has fostered new experiments and novel views on the concep- 
tual foundations of quantum mechanics, and also inspired much current research on 
coherent quantum phenomena, with quantum optical systems playing a prominent 
role. Yet, the development of quantum information had so far little impact on the 
way that quantum mechanics is taught, both at graduate and undergraduate levels. 
This tutorial is devoted to review the mathematical tools of quantum mechanics and 
to present a modern reformulation of the basic postulates which is suitable to de- 
scribe quantum systems in interaction with their environment, and with any kind of 
measuring and processing devices. 

We use Dirac notation throughout this tutorial and by system we mean a single 
given degree of freedom (spin, position, angular momentum,...) of a physical entity. 
The postulates of quantum mechanics are a list of prescriptions to summarize 

1 . how we describe the states of a physical system; 

2. how we describe the measurements performed on a physical system; 

3. how we describe the evolution of a physical system, either because of the dynamics 
or due to a measurement. 

In this section we present a picoreview of the basic postulates of quantum mechanics in 
order to introduce notation and point out both i) the implicit assumptions contained 
in the standard formulation, and ii) the need of a reformulation in terms of more 
general mathematical objects. For our purposes the postulates of quantum mechanics 
may be grouped and summarized as follows 

Postulate 1 (States of a quantum system) The possible states of a physical sys- 
tem correspond to normalized vectors li/j), (V'lV') = 1, of a Hilbert space H. Composite 
systems, either made by more than one physical object or by the different degrees of 
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freedom of the same entity, are described by tensor product Hi (g) H2 ... of the corre- 
sponding Hilbert spaces, and the overall state of the system is a vector in the global 
space. As far as the Hilbert space description of physical systems is adopted, then we 
have the superposition principle, which says that if and \ip2) are possible states 
of a system, then also any (normalized) linear combination a\tfji) + /3|V'2), a,/3 € (D, 
jap + |/3p = 1 of the two states is an admissible state of the system. 

Postulate 2 (Quantum measurements) Observable quantities are described by 
Hermitian operators X . Any hermitian operator X = X^ , admits a spectral decompo- 
sition X = xPx, in terms of its real eigenvalues x, which are the possible value of 
the observable, and of the projectors = \x){x\, Px,Px' = ^xx'Px on its eigenvectors 
X\x) = x\x), which form a basis for the Hilbert space, i.e. a complete set of orthonor- 
mal states with the properties {x\x') — Sxx' (orthonormality), and = I 

(completeness, we omitted to indicate the dimension of the Hilbert space). The prob- 
ability of obtaining the outcome x from the measurement of the observable X is given 

Px = i-e 

Px = (VIP.IV) =^(^„|P.|^)(V^|^„) =Tr[|V')(V^|P.] , (1) 

n n 

and the overall expectation value by 

{X) = {^\X\i:}=TTU){^/j\X] . 

This is the Born rule, which represents the fundamental recipe to connect the math- 
ematical description of a quantum state to the prediction of quantum theory about 
the results of an actual experiment. The state of the system after the measurement 
is the (normalized) projection of the state before the measurement on the eigenspace 
of the observed eigenvalue, i.e. 

Postulate 3 (Dynamics of a quantum system) The dynamical evolution of a 
physical system is described by unitary operators: if j'f/'o) is the state of the system 
at time to then the state of the system at time t is given by ji/'t) = V(t, to)|'0o)j with 
U{t,to)U^t,to) — {t,to)U{t,to) = I. For systems governed by a time-independent 
Hamiltonian H 7^ H{t) we have U{t,to) = exp{—iH{t — to)}. 

We will denote by L{H) the linear space of (linear) operators from H to H, which 
itself is a Hilbert space with scalar product provided by the trace operation, i.e. upon 
denoting by |^)) operators seen as elements of L{H), we have {{A\B)) — Tx[A^ B] 
(see Appendix \K\ for details on the trace operation). In the following, we are going to 
address both finite- and infinite-dimensional systems, but despite this fact we are not 
going to pay attention to the domain of operators, and to all the subtleties related to 
the geometry of infinite-dimensional Hilbert spaces. 

As it is apparent from their formulation, the postulates of quantum mechanics, as 
reported above, are about a closed isolated system. On the other hand, we are mostly 
dealing with system that interacts or have interacted with the rest of the universe, 
either during their dynamical evolution, or when subjected to a measurement. As 
a consequence, one may wonder whether a suitable modification is needed, or in 
order. This is indeed the case and the rest of his tutorial is devoted to review the 
mathematical tools of quantum mechanics and to present a modern reformulation of 
the basic postulates which is suitable to describe, design and control quantum systems 
in interaction with their environment, and with any kind of measuring and processing 
devices. 
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2 Quantum states 

2.1 Density operator and partial trace 

Suppose to have a quantum system whose preparation is not completely under con- 
trol. What we know is that the system is prepared in the state {ipk) with probability 
Pfe, i.e. that the system is described by the statistical ensemble {pk, J2kPk = 1; 

where the states {IV'fc)} are not, in general, orthogonal. The expected value of an 
observable X may be evaluated as follows 

{X) = ^Pk{X)k = ^Pk{lpk\X\ijk) = ^Pk{i^k\iPn){(Pn\X\(fp){(fp\ijk) 
k k np k 

= ^Pk{fp\lpk){i^k\'Pn){'Pn\X\Lpp) ^^{(fp\g\(pn){fn\X\ipp) 
npk np 

p 

where 

Q^'^Pk \tpk){ipk\ 

k 

is the statistical (density) operator describing the system under investigation. The 
\(Pn)'s in the above formula are a basis for the Hilbert space, and we used the trick 
of suitably inserting two resolutions of the identity I — ^„ |(p„)((^„|. The formula is 
of course trivial if the \ipk)^s are themselves a basis or a subset of a basis. 

Theorem 1 (Density operator) An operator g is the density operator associated 
to an ensemble {pk, \4'k)} is and only if it is a positive g > (hence selfadjoint) 
operator with unit trace Tr [p] = 1 . 

Proof : If g — J2kPk\^k){ipk\ is a density operator then TT[g] = J^kPk ~ ^ ^^'^ 
for any vector G H, {(p\g\ip) = J2kPk\{'P\i'k}\'^ > 0. Viceversa, if is a positive 
operator with unit trace than it can be diagonalized and the sum of eigenvalues is 
equal to one. Thus it can be naturally associated to an ensemble. □ 

As it is true for any operator, the density operator may be expressed in terms of its 
matrix elements in a given basis, i.e. g = X^np fi'nplv'n) (VpI where Qnp = ((/J„|£»|(/?p) 
is usually referred to as the density matrix of the system. Of course, the density 
matrix of a state is diagonal if we use a basis which coincides or includes the set 
of eigenvectors of the density operator, otherwise it contains off-diagonal elements. 
Different ensembles may lead to the same density operator. In this case they have the 
same expectation values for any operator and thus are physically indistinguishable. 
In other words, different ensembles leading to the same density operator are actually 
the same state, i.e. the density operator provides the natural and most fundamental 
quantum description of physical systems. 

How this reconciles with Postulate [T] dictating that physical systems are described 
by vectors in a Hilbert space? 

In order to see how it works let us first notice that, according to the postulates 
reported above, the action of "measuring nothing" should be described by the iden- 
tity operator I. Indeed the identity it is Hermitian and has the single eigenvalues 1, 
corresponding to the persistent result of measuring nothing. Besides, the eigenpro- 
jector corresponding to the eigenvalue 1 is the projector over the whole Hilbert space 
and thus we have the consistent prediction that the state after the "measurement" is 
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left unchanged. Let us now consider a situation in which a bipartite system prepared 
in the state \iPab)) G (g) is subjected to the measurement of an observable 
X — J^x^^" ^ L{Hyi), Px = \x){x\ i.e. a measurement involving only the degree 
of freedom described by the Hilbert space H^. The overall observable measured on 
the global system is thus X = X (s^lg, with spectral decomposition X = J^x^Qx^ 
= Px (^1b- The probability distribution of the outcomes is then obtained using 
the Born rule, i.e. 

Px = Tr^B [\iPab)){{^ab\Px «) Ib] . (2) 

On the other hand, since the measurement has been performed on the sole system A, 
one expects the Born rule to be valid also at the level of the single system A, and a 
question arises on the form of the object Qa which allows one to write Px = Tr^ [qa Px] 
i.e. the Born rule as a trace only over the Hilbert space Ha- Upon inspecting Eq. 
(2)) one sees that a suitable mapping \iPab)){{iPab\ — >■ 0a is provided by the partial 
trace Qa = Ttb[\'iPab)){{iPab\]- Indeed, for the operator Qa defined as the partial 
trace, we have Tr^[p^] = Ttab WiPab)) {{'^ab\] = 1 and, for any vector \(p) e Ha , 
('/'aIpaIva) = Tr4B I'ySAKv'Al <8) Ib] > 0. Being a positive, unit trace, 

operator Qa is itself a density operator according to Theorem 1. As a matter of fact, 
the partial trace is the unique operation which allows to maintain the Born rule at 
both levels, i.e. the unique operation leading to the correct description of observable 
quantities for subsystems of a composite system. Let us state this as a little theorem 

m 

Theorem 2 (Partial trace) The unique mapping |'0ab))((''A4b| Qa = /{''Pab) 
from Ha (81 Hb to Ha for which Ttab [\^ab)){{^ab\ Px <8) Ib] = Tr^ [/(V'ab) Px] is the 
partial trace f{ilJAB) ^ Qa ^ TrB [\4'ab)){{4'ab\]- 

Proof Basically the proof reduces to the fact that the set of operators on Ha is itself a 
Hilbert space L{Ha) with scalar product given by {{A\B)) = Tt[A''B]. If we consider 
a basis of operators {Mk} for L{Ha) and expand f{ipAB) — -^fcTr^[M^/(V'AB)], 
then since the map / has to preserve the Born rule, we have 

fi^AB) = ^IkTTAiMl fi^PAB)] = ^ AffcTr^B [mI ® Ib \^Ab)){{^A^ 
k k 

and the thesis follows from the fact that in a Hilbert space the decomposition on a 
basis is unique. □ 

The above result can be easily generalized to the case of a system which is initially de- 
scribed by a density operator Qab , and thus we conclude that when we focus attention 
to a subsystem of a composite larger system the unique mathematical description of 
the act of ignoring part of the degrees of freedom is provided by the partial trace. It 
remains to be proved that the partial trace of a density operator is a density operator 
too. This is a very consequence of the definition that we put in the form of another 
little theorem. 

Theorem 3 The partial traces Qa = TrB[e4B], Qb ~ Tr^[g4B] of a density operator 
Qab of a bipartite system, are themselves density operators for the reduced systems. 

Proof We have Tr^[g_4] = TrB[^?B] = Tr4B[^?AB] = 1 and, for any state \lpa) G Ha, 
Wb) G i?B, 

{'Pa\QaWa) = Tr^B [Qab \^a){Va\ «) Ib] > 

{ifB\QBWB) = Tr^B [Qab'S-A 8) \ipB){ipB\] > 0. □ 
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2.1.1 Conditional states 

From the above results it also follows that when we perform a measurement on one 
of the two subsystems, the state of the "unmeasured" subsystem after the observa- 
tion of a specific outcome may be obtained as the partial trace of the overall post 
measurement state, i.e. the projection of the state before the measurement on the 
eigenspace of the observed eigenvalue, in formula 

Qbx = — Tr^ [P^ ® la Px ® M = — TiA [Qab Px M (3) 

Vx Px 

where, in order to write the second equality, we made use of the circularity of the 
trace (see Appendix Ej) and of the fact that we are dealing with a factorized projector. 
The state Qbx will be also referred to as the "conditional state" of system B after the 
observation of the outcome x from a measurement of the observable X performed on 
the system A. 

Exercise 1 Consider a bidimensional system (say the spin state of a spin ^ particle) 
and find two ensembles corresponding to the same density operator. 

Exercise 2 Consider a spin ^ system and the ensemble {pk, IV'/cji ^ = Oi li Po = 
Pi = ^, IV'o) = |0), IV'i) = where |fc) are the eigenstates of 0-3. Write the density 
matrix in the basis made of the eigenstates of CT3 and then in the basis of cti . Then, 
do the same but for the ensemble obtained from the previous one by changing the 
probabilities to po = 3: ■ Pi = | • 

Exercise 3 Write down the partial traces of the state !')/'))= cos</) |00)) +sin0 |11)), 
where we used the notation \jk)) = \j) ® \k). 



2.2 Purity and purification of a mixed state 

As we have seen in the previous section when we observe a portion, say A, of a 
composite system described by the vector \^Pab)) G Ha®Hb, the mathematical object 
to be inserted in the Born rule in order to have the correct description of observable 
quantities is the partial trace, which individuates a density operator on Ha. Actually, 
also the converse is true, i.e. any density operator on a given Hilbert space may 
be viewed as the partial trace of a state vector on a larger Hilbert space. Let us 
prove this constructively: if p is a density operator on iJ, then it can be diagonalized 
by its eigenvectors and it can be written as q — Afc|'0/c)(V'fc|; then we introduce 
another Hilbert space K, with dimension at least equal to the number of nonzero 
eqigenvalues of g and a basis in if, and consider the vector \ip)) S H®K given 

by |(/?)) = ^/\k\il)k) ® \Qk)- Upon tracing over the Hilbert space K, we have 

Tr^[|(^))(((^|] =g VAfcAfc. \i^k){ii^k'\{ek'\eu) ^^\k\4>k){M = Q ■ 

Any vector on a larger Hilbert space which satisfies the above condition is referred to 
as a purification of the given density operator. Notice that, as it is apparent from the 
proof, there exist infinite purifications of a density operator. Overall, putting together 
this fact with the conclusions from the previous section, we are led to reformulate 
the first postulate to say that quantum states of a physical system are described by 
density operators, i.e. positive operators with unit trace on the Hilbert space of the 
system. 
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A suitable measure to quantify how far a density operator is from a projector 
is the so-called purity, which is defined as the trace of the square density operator 
fi[g] — Tr[p^] = J2k^'k' ■where the A^'s are the eigenvalues of g. Density operators 
made by a projector g = \ip){^\ have n = 1 and are referred to as pure states, 
whereas for any < 1 we have a mixed state. Purity of a state ranges in the interval 
1/d < fi < 1 where d is the dimension of the Hilbert space. The lower bound is 
found looking for the minimum of ^ ~ with the constraint Afc = 1, and 

amounts to minimize the function F — fi + 7X]fc ^k, 7 being a Lagrange multipliers. 
The solution is Xk = Vfc, i.e. the maximally mixed state g — I/d, and the 
corresponding purity is fi = I / d. 

When a system is prepared in a pure state we have the maximum possible infor- 
mation on the system according to quantum mechanics. On the other hand, for mixed 
states the degree of purity is connected with the amount of information we are miss- 
ing by looking at the system only, while ignoring the environment, i.e. the rest of the 
universe. In fact, by looking at a portion of a composite system we are ignoring the 
information encoded in the correlations between the portion under investigation and 
the rest of system: This results in a smaller amount of information about the state of 
the subsystem itself. In order to emphasize this aspect, i.e. the existence of residual 
ignorance about the system, the degree of mixedness may be quantified also by the 
Von Neumann (VN) entropy S[g] = — Tr[p \ogg] — — A„logA„, where {A„} are 
the eigenvalues of g. We have < S[g] < logd: for a pure state S'[|V')(V'|] = whereas 
— \ogd for a maximally mixed state. VN entropy is a monotone function of 
the purity, and viceversa. 

Exercise 4 Evaluate purity and VN entropy of the partial traces of the state IV')) = 
cos(/)|01)) +sin0|lO)). 

Exercise 5 Prove that for any pure bipartite state the entropies of the partial traces 
are equal, though the two density operators need not to be equal. 

Exercise 6 Take a single-qubit state with density operator expressed in terms of the 
Pauli matrices g = + riai + r2<72 + ^3(73) (Bloch sphere representation), r^ = 
Tr[fj cTfc], and prove that the Bloch vector (ri, r2, r^) should satisfies -I- -f rg < 1 
for £1 to be a density operator. 



3 Quantum measurements 

In this section we put the postulates of standard quantum measurement theory under 
closer scrutiny. We start with some formal considerations and end up with a refor- 
mulation suitable for the description of any measurement performed on a quantum 
system, including those involving external systems or a noisy environment |2l3j . 

Let us start by reviewing the postulate of standard quantum measurement theory 
in a pedantic way, i.e. by expanding Postulate [51 g denotes the state of the system 
before the measurement. 

[2.1] Any observable quantity is associated to a Hermitian operator X with spectral 
decomposition X = a;|a;)(x|. The eigenvalues are real and we assume for 
simplicity that they are nondegenerate. A measurement of X yields one of the 
eigenvalues x as possible outcomes. 

[2.2] The eigenvectors of X form a basis for the Hilbert space. The projectors Px = 
\x){x\ span the entire Hilbert space, Px = I. 

[2.3] The projectors Px are orthogonal PxPx' = ^xx'Px- It follows that P^ = Px and 
thus that the eigenvalues of any projector are and 1. 
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[2.4] (Born rule) The probability that a particular outcome is found as the measurement 
result is 

p,=Ti [P^gP,] = Ti [qP^] ^ Tr [gP,] . 

[2.5] (Reduction rule) The state after the measurement (reduction rule or projection 
postulate) is 

Qx ~ PxQPx 5 

Px 

if the outcome is x. 

[2.6] If we perform a measurement but wc do not record the results, the post-measurement 
state is given hy g = J2xP^ Sx = J2x PxQPx- 

The formulations [2.4] and [2.5] follow from the formulations for pure states, upon 
invoking the existence of a purification: 

p^ = Tr^s [Px «) Ib |V'ab))((V'ab| Px Ib] = Tr^B [\Ab)) {{-ipAsl Px Ib] 

= [QaP^] (4) 

g^^ = — Trs [Px O Ib [iIJab)) {{tlJAB[ Px ® Ib] = —Px TrB [[^^ab)) {{tlJAB[] Px 
Px Px 

= —Px Qa Px ■ (5) 

Px 

The message conveyed by these postulates is that we can only predict the spectrum of 
the possible outcomes and the probability that a given outcome is obtained. On the 
other hand, the measurement process is random, and wc cannot predict the actual 
outcome of each run. Independent! on its purity, a density operator g does not describe 
the state of a single system, but rather an ensemble of identically prepared systems. If 
we perform the same measurement on each member of the ensemble we can predict the 
possible results and the probability with which they occur but we cannot predict the 
result of individual measurement (except when the probability of a certain outcome 
is either or 1). 



3.1 Probability operator-valued measure and detection operators 

The set of postulates [2.*] may be seen as a set of recipes to generate probabilities 
and post-measurement states. We also notice that the number of possible outcomes 
is limited by the number of terms in the orthogonal resolution of identity, which itself 
cannot be larger than the dimensionality of the Hilbert space. It would however be 
often desirable to have more outcomes than the dimension of the Hilbert space while 
keeping positivity and normalization of probability distributions. In this section will 
show that this is formally possible, upon relaxing the assumptions on the mathemat- 
ical objects describing the measurement, and replacing them with more flexible ones, 
still obtaining a meaningful prescription to generate probabilities. Then, in the next 
sections we will show that there are physical processes that fit with this generalized 
description, and that actually no revision of the postulates is needed, provided that 
the degrees of freedom of the measurement apparatus are taken into account. 

The Born rule is a prescription to generate probabilities: its textbook form is the 
right term of the starred equality in [2.4]. However, the form on the left term has the 
merit to underline that in order to generate a probability it sufficient if the P^ is a 
positive operator. In fact, we do not need to require that the set of the P^^'s are pro- 
jectors, nor we need the positivity of the underlying P^ operators. So, let us consider 
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the following generalization: we introduce a set of positive operators > 0, which 
are the generalization of the and use the prescription px = TicIqIIx] to generate 
probabilities. Of course, we want to ensure that this is a true probability distribution, 
i.e. normalized, and therefore require that IIx = I, that is the positive operators 
still represent a resolution of the identity, as the set of projectors over the eigenstates 
of a selfadjoint operator. We will call a decomposition of the identity in terms of posi- 
tive operators J^x a; = I a probability operator-valued measure (POVM) and > 
the elements of the POVM. 

Let us denote the operators giving the post-measurement states (as in [2.5]) by 
Mx- We refer to them as to the detection operators. As noted above, they are no longer 
constrained to be projectors. Actually, they may be any operator with the constraint, 
imposed by [2.4] i.e. = Tr[MxgM^] = Tr[p77^]. This tells us that the POVM 
elements have the form Ux = M^M^ which, by construction, individuate a set of a 
positive operators. There is a residual freedom in designing the post-measurement 
state. In fact, since II^ is a positive operator = y/lT^ exists and satisfies the 
constraint, as well as any operator of the form = Ux \JTlx with Ux unitary. 
This is the most general form of the detection operators satisfying the constraint 
Tlx = M^Mx and corresponds to their polar decomposition. The POVM elements 
determine the absolute values leaving the freedom of choosing the unitary part. 

Overall, the detection operators Mx represent a generalization of the projectors 
Px, while the POVM elements Ux generalize F^. The postulates for quantum mea- 
surements may be reformulated as follows 

[II. 1] Observable quantities are associated to POVMs, i.e. decompositions of identity 
X^a; = I in terms of positive Ux > operators. The possible outcomes are the 
eigenvalues x of the POVM elements. 

[II. 2] The elements of a POVM are positive operators expressible as Ux = Mj Mx 
where the detection operators Mx are generic operators with the only constraint 

[II. 3] (Born rule) The probability that a particular outcome is found as the measurement 
result is = TV [M^^M]] = TV [qM^x^x] = TV [gUx]. 

[II. 4] (Reduction rule) The state after the measurement is = ^ MxqM^ if the out- 
come is X. 

[II. 5] If we perform a measurement but we do not record the results, the post-measurement 
state is given by q = J2xPxQx = Ex MxqM^. 
Since orthogonality is no longer a requirement, the number of elements of a POVM has 
no restrictions and so the number of possible outcomes from the measurement. The 
above formulation generalizes both the Born rule and the reduction rule, and says that 
any set of detection operators satisfying [II. 2] corresponds to a legitimate operations 
leading to a proper probability distribution and to a set of post-measurement states. 
This scheme is referred to as a generalized measurement. 

Of course, up to this point, this is just a formal mathematical generalization of 
the standard description of measurements given in textbook quantum mechanics, 
and few questions naturally arise: Do generalized measurements describe physically 
realizable measurements? How they can be implemented? And if this is the case, does 
it means that standard formulation is too restrictive or wrong? To all these questions 
an answer will be provided by the following sections where we state and prove the 
Naimark Theorem, and discuss few examples of measurements described by POVMs. 

3.2 The Naimark theorem 

The Naimark theorem basically says that any generalized measurement satisfying 
[II.*] may be viewed as a standard measurement defined by [2.*] in a larger Hilbert 
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space, and conversely, any standard measurement involving more than one physical 
system may be described as a generalized measurement on one of the subsystems. 
In other words, if we focus attention on a portion of a composite system where a 
standard measurement takes place, than the statistics of the outcomes and the post- 
measurement states of the subsystem may be obtained with the tools of generalized 
measurements. Overall, we have 

Theorem 4 (Naimark) For any given POVM J^x — I, > on a Hilbert 
space Ha there exists a Hilbert space H^, a state Qb = |a;i5)(a;i5| G L{Hb), a unitary 
operation U S L{Ha (S) Hg), UU^ = U^U = 1, and a projective measurement Px, 
PxP'x = Sxx'Px on Hg such that Hx — TrB[I ® Qb C/^I ® PxU]. The setup is referred 
to as a Naimark extension of the POVM. Conversely, any measurement scheme where 
the system is coupled to an ancilla and, after evolution, a projective measurement is 
performed on the ancilla may be seen as the Naimark extension of a POVM, i.e. 
one may write the Born rule Px = Tr[p^ 77a;] and the reduction rule Qa Qax = 
■^MxQaMI at the level of the system only, in terms of the POVM elements Hx = 

Ttb\^® Qb C/^I® Px U] and the detection operators Mx\fA) = {x\U\tfAi'^B))- 

Let us start with the second part of the theorem, and look at what happens when we 
couple the system under investigation to an additional system, usually referred to as 
ancilla (or apparatus), let them evolve, and then perform a projective measurement on 
the ancilla. This kind of setup is schematically depicted in Figure 1. The Hilbert space 
of the overall system is Ha ® Hb , and we assume that the system and the ancilla are 
initially independent on each other, i.e. the global initial preparation is 7? = Qa® Qb- 
We also assume that the ancilla is prepared in the pure state Qb = \uJb){^b\ since 
this is always possible, upon a suitable purification of the ancilla degrees of freedom, 
i.e. by suitably enlarging the ancilla Hilbert space. Our aim it to obtain information 
about the system by measuring an observable X on the ancilla. This is done after 
the system-ancilla interaction described by the unitary operation U. According to the 
Born rule the probability of the outcomes is given by 

Px=TrAB [UgA<S) QbU^I<»\x}{x\] = Tr^ Qa Tr^ [l<» QbU'' I(»\x){x\U] 

I " v ' 

where the set of operators 77a; = Tr^ [l I |x)(a;|J7] = {ujb\U^1(>^ PxU\uJb) is 

the object that would permit to write the Born rule at the level of the subsystem A, 
i.e. it is our candidate POVM. 




Fig. 1. Schematic diagram of a generahzed measurement. The system of interest is coupled 
to an ancilla prepared in a known state \ujb) by the unitary evolution U, and then a projective 
measurement is performed on the ancilla. 
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In order to prove this, let us define the operators G L{Ha) by their action on the 
generic vector in Ha 

M:c\lPa) = {x\U\ipA,i^B)) 

where — |(/?^) ® and the |x)'s are the orthogonal eigenvectors of X. 

Using the decomposition of — '^fe|V'fe)(V'fc| onto its eigenvectors the probability 
of the outcomes can be rewritten as 

Px^^^AB [Uqa® QbU'^^®\x){x\] =^AfcTr^B [U\^k.^^B)){{i^B,MU'^ ^®\x){x\] 

k 

= ^AfeTr^ [{x\U\iJk,^B)){{i^B,MU^\^)] -JI^feTr^ [M,\i^k){MMl] 

k k 

= Tr^ [M^QaMI] = Tr^ [qa MIM,] , (6) 

which shows that Ux — M^M^ is indeed a positive operator Va;. Besides, for any 
vector {(Pa) in Ha we have 

((p^|^AftM,|(p^) ^J2{{u:^,ipA\U^x){x\U\ipA,u:B)) 

X X 

= ((u;B:'PA\U^U\if A. OJ s)) ^1, (7) 

and since this is true for any \lpa) we have J^x ^l^x = I- Putting together Eqs. ^ 
and (O we have that the set of operators H^ = M^M^ is a POVM, with detection 
operators M^- In turn, the conditional state of the system A, after having observed 
the outcome x, is given by 

Qax = — Tre [U QA®\i^B){<-JB\U'^ ^® Px] ^ — ^Xk{x\U\i;k,i^B)){{^B.MU'^\x) 

Px Px , 

k 

= —M^QaMI (8) 

Px 

This is the half of the Naimark theorem: if we couple our system to an ancilla, let 
them evolve and perform the measurement of an observable on the ancilla, which 
projects the ancilla on a basis in iJ^, then this procedure also modify the system. 
The transformation needs not to be a projection. Rather, it is adequately described 
by a set of detection operators which realizes a POVM on the system Hilbert space. 
Overall, the meaning of the above proof is twofold: on the one hand we have shown 
that there exists realistic measurement schemes which are described by POVMs when 
we look at the system only. At the same time, we have shown that the partial trace of 
a spectral measure is a POVM, which itself depends on the projective measurement 
performed on the ancilla, and on its initial preparation. Finally, we notice that the 
scheme of Figure 1 provides a general model for any kind of detector with internal 
degrees of freedom. 

Let us now address the converse problem: given a set of detection operators Mx 
which realizes a POVM M^M^ — I, is this the system-only description of an 
indirect measurement performed a larger Hilbert space? In other words, there exists 
a Hilbert space Hg, a state Qb = |a;i5)(a;B| G L{Hb), a unitary U G L{Ha Hg), 
and a projective measurement = |a;)(a::| in Hg such that M^lfA) — {x\U\ipa,'-^b)) 
holds for any \(Pa) & Ha and H^ — (wslC/^'I PxU\ujb)'^ The answer is positive and 
we will provide a constructive proof. Let us take H^ with dimension equal to the 
number of detection operators and of POVM elements, and choose a basis |a;) for Hg, 
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which in turn individuates a projective measurement. Then we choose an arbitrary 
state \lOb) € Hb and define the action of an operator U as 

U\iPa) O I^b) = ^Mj; \ipa) O |a;) 

X 

where \ipA € Ha is arbitrary. The operator U preserves the scalar product 
{{uJs^^'AW'^m^A^uJs)) = ^{^'a\mImM{x'\x) = ^{ip'^lMlMM = ivM 

xx' X 

and so it is unitary in the one-dimensional subspace spanned by \lob)- Besides, it 
may be extended to a full unitary operator in the global Hilbert space Ha (81 H^ , eg 
it can be the identity operator in the subspace orthogonal to \oJb)- Finally, for any 
\(Pa) € Ha, we have 

{x\U\(Pa,cob)) ='^Mro'\ipA){x\x') = Mx\ipa) , 

x' 

and 

{ipA\n^\ipA) = {ipA\M^^M^\ipA) = {{LOB,ipA\U''I(^ P^U\ipA,iOB)) , 

tha.t is, Hx = {lJb\U^'^® PxU\uib). □ 
This completes the proof of the Naimark theorem, which asserts that there is a 
one-to-one correspondence between POVM and indirect measurements of the type 
describe above. In other words, an indirect measurement may be seen as the physical 
implementation of a POVM and any POVM may be realized by an indirect measure- 
ment. 

The emerging picture is thus the following: In measuring a quantity of interest on 

a physical system one generally deals with a larger system that involves additional 
degrees of freedom, besides those of the system itself. These additional physical en- 
tities are globally referred to as the apparatus or the ancilla. As a matter of fact, 
the measured quantity may be always described by a standard observable, however 
on a larger Hilbert space describing both the system and the apparatus. When we 
trace out the degrees of freedom of the apparatus we are generally left with a POVM 
rather than a PVM. Conversely, any conceivable POVM, i.e. a set of positive oper- 
ators providing a resolution of identity, describe a generalized measurement, which 
may be always implemented as a standard measurement in a larger Hilbert space. 
Before ending this Section, few remarks are in order; 

Rl The possible Naimark extensions are actually infinite, corresponding to the in- 
tuitive idea that there are infinite ways, with an arbitrary number of ancillary 
systems, of measuring a given quantity. The construction reported above is some- 
times referred to as the canonical extension of a POVM. The Naimark theorem 
just says that an implementation in terms of an ancilla-based indirect measure- 
ment is always possible, but of course the actual implementation may be diflferent 
from the canonical one. 

R2 The projection postulate described at the beginning of this section, the scheme of 
indirect measurement, and the canonical extension of a POVM have in common 
the assumption that a nondemolitive detection scheme takes place, in which the 
system after the measurement has been modified, but still exists. This is sometimes 
referred to as a measurement of the first kind in textbook quantum mechanics. 
Conversely, in a demolitive measurement or measurement of the second kind, the 
system is destroyed during the measurement and it makes no sense of speaking of 
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the state of the system after the measurement. Notice, however, that for demoli- 
tive measurements on a field the formalism of generalized measurements provides 
the framework for the correct description of the state evolution. As for example, 
let us consider the detection of photons on a single-mode of the radiation field. 
A demolitve photodetector (as those based on the absorption of light) realizes, in 
ideal condition, the measurement of the number operator a' a without leaving any 
photon in the mode . If g = J^np Qnp\n){p\ is the state of the single-mode radia- 
tion field a photodetector of this kind gives a natural number n as output, with 
probability p„ — Qnn, whereas the post-measurement state is the vacuum |0)(0| 
independently on the outcome of the measurement. This kind of measurement is 
described by the orthogonal POVM i7„ — \n){n\, made by the eigenvectors of the 
number operator, and by the detection operator A/„ = |0)(rt|. The proof is left as 
an exercise. 

R3 We have formulated and proved the Naimark theorem in a restricted form, suitable 
for our purposes. It should be noticed that it holds in more general terms, as for 
example with extension of the Hilbert space given by direct sum rather than tensor 
product, and also relaxing the hypothesis [4]. 

3.2.1 Conditional states in generalized measurements 

If we have a composite system and we perform a projective measurement on, 
say, subsystem A, the conditional state of the unmeasured subsystem B after 
the observation of the outcome x is given by Eq. i.e. it is the partial trace 
of the projection of the state before the measurement on the eigenspace of the 
observed eigenvalue. One may wonder whether a similar results holds also when 
the measurement performed on the subsystem a A is described by a POVM. The 
answer is positive and the proof may be given in two ways. The first is based on 
the observation that, thanks to the existence of a canonical Naimark extension, 
we may write the state of the global system after the measurement as 

0ABX = —M.^ (g) Ig Qab (g) Ib , 

Px 

and thus the conditional state of subsystem B is the partial trace Qbx — Tr^ [gAfsa;] 
i.e. 

Qbx = —T^AW^ms qab Mlms] = — Tr^[(?AB MlM^ms] = —^^a[qab n^mu] , 

Px Px Px 

where again we used the circularity of partial trace in the presence of factorized 
operators. A second proof may be offered invoking the Naimark theorem only to 
ensure the existence of an extension, i.e. a projective measurement on a larger 
Hilbert space He <E) Hji, which reduces to the POVM after tracing over He- In 
formula, assuming that Px G L{Hc <E) Ha) is a projector and a E L{Hc) a density 
operator 

Qbx = — TrcA [Px ® Ib Qab ®(jPx®'^b\ = —T^^CA [Qab ® (J Px ® Is] 
Px Px 

= — Tr^ [qabHx ® Ib] • 

Px 

3.3 Joint measurement of non commuting observables 

A common statement about quantum measurements says that it is not possible to 
perform a joint measurement of two observables Qa and Pa of a given system A 
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if they do not commute, i.e. [Qa^Pa] 7^ 0. This is related to the impossibility of 
finding any common set of projectors on the Hilbert space Ha of the system and to 
define a joint observable. On the other hand, a question arises on whether common 
projectors may be found in a larger Hilbert space, i.e. whether one may implement a 
joint measurement in the form of a generalized measurement. The answer is indeed 
positive |5I6| : This Section is devoted to describe the canonical implementation of joint 
measurements for pair of observables having a (nonzero) commutator [Qa, Pa] = cI 7^ 
proportional to the identity operator. 

The basic idea is to look for a pair of commuting observables [JCis, I^b] = in a 
larger Hilbert space ® Hg which trace the observables Pa and Qa, i-C. which have 
the same expectation values 

(Xab) = Ttab[Xab Qa 8) Qb] = Tr^[QA Qa] = (Qa) 
(Yab) = TVab[Yab Qa ® Qb] = Tr^[P^ Qa] = {Pa) (9) 

for any state Qa G Ha of the system under investigation, and a fixed suitable prepa- 
ration Qb G Hb of the system B. A pair of such observables may be found upon 
choosing a replica system B, identical to A, and considering the operators 

Xab = Qa®^b+^a®Qb 

Yab = Pa®^b-Ia®Pb (10) 

where Qb and Pg are the analogue of Qa and Pa for system B. The operators in Eq. 
(fTU)l . taken together a state Qb & Hb satisfying 

Tr^ [Qb Qb] = Tr^ [Pb Qb] = , (11) 

fulfill the conditions in Eq. i.e. realize a joint generalized measurement of the 
noncommuting observables Qa and Pa- The operators Xab and Yab are Hermitian by 
construction. Their commutator is given by 

[Xab, Yab] = [Qa, Pa] <S) Ib - ^a [Qb, Pb] ^ ■ (12) 

Notice that the last equality, i.e. the fact that the two operators commute, is valid 
only if the commutator [Qa,Pa] = cI is proportional to the identity. More general 
constructions are needed if this condition does not hold [7; . 

Since the [JCibj^b] = the complex operator Zab = Xab + is normal i.e. 
[Zab, zIb] — 0. For normal operators the spectral theorem holds, and we may write 

Zab=J2^P^' = k»((^l Zab[z)) = z\z)) (13) 

z 

where z e C, and P^ are orthogonal projectors on the eigenstates |z)) = [zJiab of Zab- 
The set \Pz\ represents the common projectors individuating the joint observable 
Zab- Each run of the measurement returns a complex number, whose real and imag- 
inary parts correspond to a sample of the Xab and Yab values, aiming at sampling 
Qa and Pa- The statistics of the measurement is given by 

Vz{z) = TyabXqa ® Qb Pz] = Tr^^A H,] (14) 
where the POVM H^ is given by 



H.^TTBilA^QBPz 



(15) 
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The mean values {Xab) = (Qa) and (Yab) — (Pa) are the correct ones by construction, 
where by saying " correct" we intend the mean values that one would have recorded by 
measuring the two observables Qa and Pa separately in a standard (single) projective 
measurement on Qa- On the other hand, the two marginal distributions 

Px{x) = J dypz{x + iy) p^iy) = Jdxpz{x + iy) , 

need not to reproduce the distributions obtained in single measurements. In particu- 
lar, for the measured variances {AXa^) ~ {^Ib) ~ {Xab)"^ and {AYab) one obtains 

{AXIb) = Tr [(Ql ®Js+^a®QI + '2Qa® Qb) Qa ® Qb] - (Qa)^ 

{AYaD ^ {API) + {Pb)' (16) 

where we have already taken into account that (Qb) — (Pb) — 0. As it is apparent 
from Eqs. ([TB| the variances of Xab and Yab are larger than those of the original, non 
commuting, observables Qa and Pa- 

Overall, we may summarize the emerging picture as follows: a joint measurement 
of a pair of non commuting observables corresponds to a generalized measurement 
and may be implemented as the measurement of a pair of commuting observables on 
an enlarged Hilbert space. Mean values are preserved whereas the non commuting 
nature of the original observables manifests itself in the broadening of the marginal 
distributions, i.e. as an additional noise term appears to both the variances. The 
uncertainty product may be written as 



{AX!b){AY,1) = {AQl){APj) + {AQl){Pl) + {Ql){APl) + {Ql){Pl) , 



> \\[Qa.Pa]\^ + {AQl){Pl) + {Ql){APl) + {Ql){Pl) , (17) 

where the last three terms are usually referred to as the added noise due to the joint 
measurement. If we perform a joint measurement on a minimum uncertainty state 
(MUS, see Appendix [B|) for a given pair of observables (e.g. a coherent state in the 
joint measurement of a pair of conjugated quadratures of the radiation field) and use 
a MUS also for the preparation of the rephca system (e.g. the vacuum), then Eq. p7| 
rewrites as 

{AXI,){AY^,) = \[Qa,Pa]\^. (18) 

This is four times the minimum attainable uncertainty product in the case of a mea- 
surement of a single observable (see Appendix IbI) . In terms of rms' AX = y^{AX^ 
we have a factor 2, which is usually referred to as the 3 dB of added noise in joint 
measurements. The experimental realization of joint measurements of non commuting 
observables has been carried out for conjugated quadratures of the radiation field in 
a wide range of frequencies ranging from radiowaves to the optical domain, see e.g. 
0. 



3.4 About the so-called Heisenberg principle 

Let us start by quoting Wikipedia about the Heisenberg principle [S] 

Published by Werner Heisenberg in 1927, the principle implies that it is impos- 
sible to simultaneously both measure the present position while "determining" 
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the future momentum of an electron or any other particle with an arbitrary 
degree of accuracy and certainty. This is not a statement about researchers' 
ability to measure one quantity while determining the other quantity. Rather, 
it is a statement about the laws of physics. That is, a system cannot be de- 
fined to simultaneously measure one value while determining the future value 
of these pairs of quantities. The principle states that a minimum exists for the 
product of the uncertainties in these properties that is equal to or greater than 
one half of the reduced Planck constant. 

As is it apparent from the above formulation, the principle is about the precision 
achievable in the measurement of an observable and the disturbance introduced by 
the same measurement on the state under investigation, which, in turn, limits the 
precision of a subsequent measurement of the conjugated observable. The principle, 
which has been quite useful in the historical development of quantum mechanics, has 
been inferred from the analysis of the celebrated Heisenberg' gedanken experiments, 
and thus is heuristic in nature. However, since its mathematical formulation is related 
to that of the uncertainty relations (see Appendix [B|) . it is often though as a theorem 
following from the axiomatic structure of quantum mechanics. This is not the case: 
here we exploit the formalism of generalized measurements to provide an explicit 
example of a measurement scheme providing the maximum information about a given 
observable, i.e. the statistics of the corresponding PVM, while leaving the state under 
investigation in an eigenstate of the conjugated observable. 

Let us consider the two noncommuting observables [A, B] = clL and the set of 
detection operators Ma = \b){a\ where \a) and |6) are eigenstates of A and B respec- 
tively, i.e. A\a) = a\a), B\b) = b\b). According to the Naimark theorem the set of 
operators {Ma} describe a generalized measurement (e.g. an indirect measurement 
as the one depicted in Fig. 1) with statistics Pa — Tr[gi ila] described by the POVM 
77a = M^Ma = \a){a\ and where the conditional states after the measurement are 
given by ga = ^MagM^ = \b){b\. In other words, the generalized measurement 
described by the set {Ma} has the same statistics of a Von-Neumann projective mea- 
surement of the observable A, and leave the system under investigating in an eigen- 
state of the observable B, thus determining its future value with an arbitrary degree 
of accuracy and certainty and contrasting the formulation of the so-called Heisenberg 
principle reported above. An explicit unitary realization of this kind of measurement 
for the case of position, as well as a detailed discussion on the exact meaning of 
the Heisenberg principle, and the tradeoff between precision and disturbance in a 
quantum measurement, may be found in [10) . 



3.5 The quantum roulette 

Let us consider K projective measurements corresponding to K nondegenerate isospec- 
tral observables Xk, k = 1,...,K in a Hilbert space H, and consider the following 
experiment. The system is sent to a detector which at random, with probability 
Zfc, Zk = 1, perform the measurement of the observable Xk- This is known as the 
quantum roulette since the observable to be measured is chosen at random, eg accord- 
ing to the outcome of a random generator like a roulette. The probability of getting 
the outcome x from the measurement of the observable X^ on a state g G L{H) is 

given by p^x^ = Tr[gi Pi''^], pj*^^ = \x)kk{x\, and the overall probability of getting the 
outcome x from our experiment is given by 

P^ = T. '^f^P^'^ - E ^fcTr[^. Pi^-)] = Tr[^, ^ z^Pi'^] = Tr[f,77,] , 
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where the POVM describing the measurement is given by IIx — ^kPx ■ This is 
indeed a POVM and not a projective measurement since 

kk' 

Again, we have a practical situation where POVMs naturally arise in order to describe 
the statistics of the measurement in terms of the Born rule and the system density 
operator. A Naimark extension for the quantum roulette may be obtained as follows. 
Let us consider an additional probe system described by the Hilbert space Hp of 
dimension K equal to the number of measured observables in the roulette, and the 
set of projectors Qx = J2k ^^^'^ ® \0k){0k\ where {16*;:)} is a basis for Hp. Then, upon 
preparing the probe system in the superposition \ujp) = ^/zk\0k) we have that 

Px = TtspIq® \i^p){<^p \ Qx] and, in turn, H^ = Trp[Is ® |wp)(wp| Q^] = T.k^kPx^'' ■ 
The state of the system after the measurement may be obtained as the partial trace 

Qx = — Trp [QxQ® \ujp){LOp\ Qx] 

Px 

Pi'^ ® \0k){0k\g<E> \ujp){ujp\ pp <E) \ek'){ek'\ 

''■^ k k' 

= -Y.^kPi''^p'^''- 

Notice that the presented Naimark extension is not the canonical one. 

Exercise 7 Prove that the operators Qx introduced for the Naimark extension of 
the quantum roulette, are indeed projectors. 

Exercise 8 Take a system made by a single qubit system and construct the canonical 
Naimark extension for the quantum roulette obtained by measuring the observables 
a a = cosacTi + sinQ!a2, where ai and 1J2 are Pauli matrices and a £ [0, tt] is chosen 
at random with probability p{a) = tt^^. 



4 Quantum operations 

In this Section we address the dynamical evolution of quantum systems to see whether 
the standard formulation in terms of unitary evolutions needs a suitable generaliza- 
tion. This is indeed the case: we will introduce a generalized description and see how 
this reconciles with what we call Postulate [3] in the Introduction. We will proceed 
in close analogy with what we have done for states and measurements. We start by 
closely inspecting the physical motivations behind any mathematical description of 
quantum evolution, and look for physically motivated conditions that a map, intended 
to transform a quantum state into a quantum state, from now on a quantum operation, 
should satisfy to be admissible. This will lead us to the concept of complete positivity, 
which suitably generalizes the motivations behind unitarity. We then prove that any 
quantum operation may be seen as the partial trace of a unitary evolution in a larger 
Hilbert space, and illustrate a convenient form, the so-called Kraus or operator-sum 
representation, to express the action of a quantum operation on quantum states. 

By quantum operation we mean a map Q ^ £{q) transforming a quantum state g 
into another quantum state £{q). The basic requirements on £ to describe a physically 
admissible operations are those captured by the request of unitarity in the standard 
formulation, i.e. 
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Ql The map is positive and trace-preserving, i.e. £{g)] > (hence selfadjoint) and 
Tr[£(p)] = Tr[p] = 1. The last assumption may be relaxed to that of being trace 
non-increasing < Tr[£((3)] < 1 in order to include evolution induced by mea- 
surements (see below). 

Q2 The map is linear £{J2kPk9k) = J^kP^^i^k), i-e. the state obtained by applying 
the map to the ensemble {pk,(>k} is the ensemble {pk,£(gk)}- 

Q3 The map is completely positive (CP), i.e. besides being positive it is such that if we 
introduce an additional system, any map of the form £ 1 acting on the extended 
Hilbert space is also positive. In other words, we ask that the map is physically 
meaningful also when acting on a portion of a larger, composite, system. As we 
will see, this request is not trivial at all, i.e. there exist maps that are positive but 
not completely positive. 

4.1 The operator-sum representation 

This section is devoted to state and prove a theorem showing that a map is a quantum 
operation if and only if it is the partial trace of a unitary evolution in a larger Hilbert 
space, and provides a convenient form, the so-called Kraus decomposition or operator- 
sum representation |11I12| . to express its action on quantum states. 

Theorem 5 (Kraus) A map £ is a quantum operation i.e. it satisfies the require- 
ments Q1-Q3 if and only if is the partial trace of a unitary evolution on a larger 
Hilbert space with factorized initial condition or, equivalently, it possesses a Kraus 
decomposition i. e. its action may be represented as £{g) = MkgMj, where {Mk} 
is a set of operators satisfying M^Mk — I. 

Proof The first part of the theorem consists in assuming that £{g) is the partial trace 
of a unitary operation in a larger Hilbert space and prove that it has a Kraus decom- 
position and, in turn, it satisfies the requirements Q1-Q3. Let us consider a physical 
system A prepared in the quantum state and another system B prepared in the 
state gg. A and B interact through the unitary operation U and we are interested 
in describing the effect of this interaction on the system A only, i.e. we are looking 
for the expression of the mapping g'^ — £{gA) induced by the interaction. This 

may be obtained by performing the partial trace over the system B of the global AB 
system after the interaction, in formula 

f(^^^) =Tr« [Ug^^gsU^] = ^p,Tr^ [C/ ® |0,)(0,|C/t] 

S 

St k 

where we have introduced the operator Mk = ^/ph{ft\U\9s), with the polyindex k = st 
obtained by a suitable ordering, and used the spectral decomposition of the density 
operator g^ = '^sPs\0s){0s\. Actually, we could have also assumed the additional 
system in a pure state |ws), since this is always possible upon invoking a purification, 
i.e. by suitably enlarging the Hilbert space. In this case the elements in the Kraus 
decomposition of our map would have be written as {ipt\U\ujB). The set of operators 
{Mk} satisfies the relation 

k st s 

Notice that the assumption of a factorized initial state is crucial to prove the existence 
of a Kraus decomposition and, in turn, the complete positivity. In fact, the dynamics 
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£{Qa) = Tr^ [[/ Qab U'^\ resulting from assuming an initially correlated preparation 
Qab needs not be so. In this case, the dynamics can properly be defined only on a 
subset of initial states of the system. Of course, the map can be extended to all possible 
initial states is by linearity, but the extension may not be physically realizable, i.e. 
may be not completely positive or even positive [T5] . 

We now proceed to show that for map of the form (Kraus decomposition) the 
properties Ql-QS hold. Preservation of trace and of the Hermitian character, as well 
as linearity, are guaranteed by the very form of the map. Positivity is also ensured, 
since for any positive operator Oa G L{Ha) and any vector G Ha we have 

{ipA\S{OA)WA) = {ipA\Y,MkOAMl\ipA) = {ipA[TrB[UOA®QB U^]Wa) 

k 

= Tvab[U^Wa){Va\®WOa®Qb]>(^ ^Oa^ QB-yWA) ■ 

Therefore it remains to be proved that the map is completely positive. To this aim 
let us consider a positive operator Oac S L{Ha <S) He) and a generic state |Vkc)) on 
the same enlarged space, and define 

\oJk)) = -i=Mk <E) Icli^Ac)) Nk = {{^JAc\MlMk (E) IcIV-Ac)) > . 

V ^^fc 

Since Oac is positive we have 

{{lpAc\{Ml(g)Ic)OAciMk(E>IcMAc)) = Nk{{L0k\OAcH)) >0 

and therefore ((i/'acI^ Ic(O^c)IV'ac)) = J2k ^k{{i^k\OAc\^^k)) > 0, which proves 
that for any positive Oac also £ (8) 1c{Oac) is positive for any choice of He, i.e. £ is 
a CP-map. 

Let us now prove the second part of the theorem, i.e. we consider a map £ : 
L{Ha) — > L{Ha) satisfying the requirements Ql-QS and show that it may be written 
in the Kraus form and, in turn, that its action may be obtained as the partial trace 
of a unitary evolution in a larger Hilbert. We start by considering the state \<f)) = 

Sfc \0k)® \dk) G Ha® Ha and define the operator Qaa — £ {{'^\) ■ From the 
complete positivity and trace preserving properties of £ we have that Tr[£i^^] — 1, 
and Qaa ^ 0, i.e. Qaa is a density operator. Besides, this establishes a one-to-ne 
correspondence between maps L{Ha) L{Ha) and density operators in L{Ha) €5 
L{Ha) which may be proved as follows: for any — X^fcV'fcl^fc) S Ha define j-^) = 
SfcV'fcl^fc) and notice that 

{4'\0aaW = ^{^p\J2£m){eiMek){ei\ |^) = v,>fef(|0fe>(e/l) = 

kl kl 

where we used linearity to obtain the last equality. Then define the operators Mk\ip) = 
\/c?pfe(V'|wfe)), where \u!k)) are the eigenvectors of Qaa = J2kPk\'-^k)){{^k\- this is a 
linear operator on Ha and we have 

k k 

for all pure states. Using again linearity we have that £{g) = J2k ^^^qMI also for 
any mixed state. It remains to be proved that a unitary extension exists, i.e. to 
prove that for any map on L{Ha) which satisfies Ql-QS, and thus possesses a Kraus 
decomposition, there exist: i) a Hilbert space Hb, ii) a state \u)b) G Hb, iii) a unitary 



20 



Will be inserted by the editor 



U e L(-H"a O Hb) such that Sig^) = Trs[U Qa ® \uj b) {oj bW^] for any G L{Ha). To 
this aim we proceed as we did for the proof of the Naimark theorem, i.e. we take an 
arbitrary state \uJb) & Hg, and define an operator U trough its action on the generic 

(Pa) ® \ujb) e -ff a ® Hb, U\lpa) ® \ojb) = J2k Wa) ® \6k), where the |^fc)'s are a 
basis for Hb- The operator U preserves the scalar product 

{{uJb,<p'AUM'Pa,uJb)) = J2{v'A\Ml,Mk\ipA){ek'm = ^(v'JMtMfe|(^,) = {^'M 

kk' k 

and so it is unitary in the one-dimensional subspace spanned by |wb). Besides, it may 
be extended to a full unitary operator in the global Hilbert space H^ (S) Hb , eg it can 
be the identity operator in the subspace orthogonal to \oJb)- Then, for any in H^ 
we have 

TVb [C/^^O |a;B)(wB I [/■•■] =^p«TrB [U\ips){i^s\® \i^B){i^B\U''] 

S 

= Y^PsTvb [Mk\i;,){i;s\Ml 

skk' 

= J2p^ Mk\^s){tps\ Ml = Mkg^Ml □ 

sk k 

The Kraus decomposition of a quantum operation generalizes the unitary description 
of quantum evolution. Unitary maps are, of course, included and correspond to maps 

whose Kraus decomposition contains a single elements. The set of quantum operations 
constitutes a semigroup, i.e. the composition of two quantum operations is still a 
quantum operation: 

kl fcl/C2 fe 

where we have introduced the polyindex k. Normalization is easily proved, since 

Y,k MlMk = Efeife2 ^ia^^^^i^^^^fei^^fea^ = I- On other hand, the existence of 
inverse is not guaranteed: actually only unitary operations are invertible. 

The Kraus theorem also allows us to have a unified picture of quantum evohition, 
either due to an interaction or to a measurement. In fact, the modification of the state 
in the both processes is described by a set of operators satisfying J2k M^Mk = I. 
In this framework, the Kraus operators of a measurement are what we have referred 
to as the detection operators of a POVM. 



4.1.1 The dual map and the unitary equivalence 

Upon writing the generic expectation value for the evolved state £{g) and exploiting 
both linearity and circularity of trace we have 

{X) = Tr[£{g)X] = Y^TiiMkgMlx] = ^ Tr[eM^^XMfe] = Ti[g£^ {X)\ , 

k k 

where we have defined the dual map E'^ {X) ~ J^k^lXMk which represents the 
"Heisenberg picture" for quantum operations. Notice also that the elements of the 
Kraus decomposition = {ipk\U\LOB) depend on the choice of the basis used to 
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perform the partial trace. Change of basis cannot have a physical effect and this 
means that the set of operators 

TVfc = (0fc|i7|w«) - ^(0fe|(^3)((^,|C/K) = ^14.M, , 

s s 

where the unitary V G L{Hb) describes the change of basis, and the original set Mk 
actually describe the same quantum operations, i.e. Yl,k^^S^k ~ ^k-^kQMl, Vg. 
The same can be easily proved for the system B prepared in mixed state. The origin 
of this degree of freedom stays in the fact that if the unitary U on (g) Hg and 
the state \uJb) G Hb realize an extension for the map £ : L{Ha) — >■ L{Ha) then any 
unitary of the form (1(g) y)J7 is a unitary extension too, with the same ancilla state. A 
quantum operation is thus identified by an equivalence class of Kraus decompositions. 
An interesting corollary is that any quantum operation on a given Hilbert space of 
dimension d may be generated by a Kraus decomposition containing at most d? ele- 
ments, i.e. given a Kraus decomposition £{q) = M^gM^ with an arbitrary number 
of elements, one may exploit the unitary equivalence and find another representation 
£{g) = J2k -^^qNI with at most elements. 



4.2 The random unitary map and the depolarizing channel 

A simple example of quantum operation is the random unitary map, defined by the 
Kraus decomposition £{q) = ^^PkUkQUj., i.e. Mk = ^/PkUk and UkUl = 1. This 
map may be seen as the evolution resulting from the interaction of our system with 
another system of dimension equal to the number of elements in the Kraus decompo- 
sition of the map via the unitary V defined by V^IV'a) (X" I^^b) = J2k y^Uk\ipA) (Xi 
\9k) being a basis for Hg which includes \u!b)- If "we do not look" at the system 
B and trace out its degree of freedom the evolution of system A is governed by the 
random unitary map introduced above. 

Exercise 9 Prove explicitly the unitarity of V. 

The operator-sum representation of quantum evolutions have been introduced, and 
finds its natural application, for the description of propagation in noisy channels, i.e. 
the evolution resulting from the interaction of the system of interest with an external 
environment, which generally introduces noise in the system degrading its coherence. 
As for example, let us consider a qubit system (say, the polarization of a photon), on 
which we have encoded binary information according to a suitable coding procedure, 
traveling from a sender to a receiver. The propagation needs a physical support (say, 
an optical fiber) and this unavoidably leads to consider possible perturbations to 
our qubit, due to the interaction with the environment. The resulting open system 
dynamics is usually governed by a Master equation, i.e. the equation obtained by 
partially tracing the Schroedinger (Von Neumann) equation governing the dynamics 
of the global system, and the solution is expressed in form of a CP-map. For a qubit 
Q in a noisy environment a quite general description of the detrimental effects of the 
environment is the so-called depolarizing channel [T] , which is described by the Kraus 
operator Mq = ^/T^^ao, Mk = (7k, k = 1,2, 3, i.e. 

f(g) = (l-7)e+^^afe^.afe 0<7<1. 
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The depolarizing channel may be seen as the evolution of the qubit due to the inter- 
action with a four-dimensional system through the unitary 

^ fe=l 

\9k) being a basis which includes |we)- From the practical point view, the map de- 
scribes a situation in which, independently on the underlying physical mechanism, we 
have a probability 7/8 that a perturbation described by a Pauli matrix is applied to 
the qubit. If we apply ai we have the so-called spin-flip i.e. the exchange |0) -H- 
whereas if we apply 0-3 we have the phase-flip, and for a2 we have a specific combi- 
nation of the two effects. Since for any state of a qubit Q + J2k ^kQO'k = 21 the action 
of the depolarizing channel may be written as 

£{g) = {l-j)g+l{2I-g) = - ^i)q = PQ + {1 - p)^ , 

where p= 1 — I7, i.e. — ^ < p < 1. In other words, we have that the original state g is 
sent to a linear combination of itself and the maximally mixed state | , also referred 

to as the depolarized state. 

Exercise 10 Express the generic qubit state in Bloch representation and explicitly 
write the effect of the depolarizing channel on the Bloch vector. 

Exercise 11 Show that the purity of a qubit cannot increase under the action of the 
depolarizing channel. 



4.3 Transposition and partial transposition 

The transpose T{X) = X'^ of an operator X is the conjugate of its adjoint X"^ = 
[X'^Y = {X*Y . Upon the choice of a basis we have X = ^nk-^nk\Sn){Sk\ and 
thus X'^ = J2nk^nk\dk){0n\ = J2nk ^kn\On) {Ok\- Transposition docs not change the 
trace of an operator, neither its eigenvalues. Thus it transform density operators 
into density operators: Tr[g] = Tr[^'^] = 1 > if ^ > 0. As a positive, trace 
preserving, map it is a candidate to be a quantum operation. On the other hand, 
we will show by a counterexample that it fails to be completely positive and thus it 
does not correspond to physically admissible quantum operation. Let us consider a 
bipartite system formed by two qubits prepared in the state \(p)) — |00)) + 

We denote hy g'^ = KS) T{g) the partial transpose of g i.e. the operator obtained by 
the application of the transposition map to one of the two qubits. We have 



/I r 
0000 
0000 

Vi 1, 



i(|0)(0| ® |0)(0| + |1)(1| ® + |0)(1| ® |0)(1| + |1)(0| ® |1)(0|) 
i(|0)(0| ® |0)(0| + |1)(1| ® |1)(1| + |0)(1| ® |1)(0| + |1>(0| ® |0>(1 



/I 0^ 
10 
10 

Vo 1, 
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Using the last expression it is straightforward to evaluate the eigenvalues of , which 
are +5 (multiplicity three) and — ^. In other words I ® T is not a positive map and 
the transposition is not completely positive. Notice that for a factorized state of the 
form = Qa® Qb wc have I (g) T{qab) = Qa ^ Q% > i.e. partial transposition 

preserves positivity in this case . 

Exercise 12 Prove that transposition is not a CP-map by its action on any state 
of the form \ip)) = ^J2k \V'k) \^k)- Hint: the operator I (g) T{\ip)){{ip\) = E is the 
so-called swap operator since it "exchanges" states as £'(|V')a \v)b) = \v)a <S) |V')b- 

5 Conclusions 

In this tutorial, we have addressed the postulates of quantum mechanics about states, 
measurements and operations. We have reviewed their modern formulation and intro- 
duced the basic mathematical tools: density operators, POVMs, detection operators 
and CP-maps. We have shown how they provide a suitable framework to describe 
quantum systems in interaction with their environment, and with any kind of mea- 
suring and processing devices. The connection with the standard formulation have 
been investigated in details building upon the concept of purification and the Theo- 
rems of Naimark and Stinespring/Kraus-Choi-Sudarshan. 

The framework and the tools illustrated in this tutorial are suitable for the pur- 
poses of quantum information science and technology, a field which has fostered new 
experiments and novel views on the conceptual foundation of quantum mechanicss, 
but has so far little impact on the way that it is taught. We hope to contribute in 
disseminating these notions to a larger audience, in the belief that they are useful for 
several other fields, from condensed matter physics to quantum biology. 
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A Trace and partial trace 

The trace of an operator O is a scalar quantity equal to sum of diagonal elements in a 
given basis Ti[0] = J2n(fn\^\'^"-) ■ "^^^ trace is invariant under any change of basis, 
as it is proved by the following chain of equalities 

Y,{0n\O\9,,) = Y,{0n\Vk){V^k\O\^j){Vj\en) = Y,(Vj\On) {9,M) {^k\0\^j) 
n njk njk 

= '^{fj\^k){'Pk\0\ip.j) = J^ifklOlfk) , 
jk k 

where we have suitably inserted and removed resolutions of the identity in terms of 
both basis {|^n)} and {| (/?„)}. As a consequence, using the basis of eigenvectors of 
O, Tr[0] = X^n*-*"' being the eigenvalues of O. Trace is a linear operation, i.e. 
Tr[Oi + O2] = Tr[Oi] + Tr[02] and Tr[AO] = ATr[0] and thus 9Tr[0] = Ti[dO] for 
any derivation. The trace of any "ket-bra" Tr[|V'i)('02|] is obtained by "closing the 
sandwich" Tr[|'0i)('02|] = (V'2|'0i)i i^i fact upon expanding the two vectors in the same 
basis and taking the trace in that basis Tr[|?/'i)(V'2|] = Snfci V'ifcV'2i(^n|^fc)(^i|^n) = 
V'inV'2n — (V'2|'*/'i)- Other properties are summarized by the following theorem. 

Theorem 6 For the trace operation the following properties hold 

i) Given any pair of operators Tt[AiA2] = Tt[A2Ai] 
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ii) Given any set of operators Ai, we Tr[^i^2^3---^jv] = Tr[j42A3...AjvAi] = 
Tr[yl3yl4...Ai A2] = ... (circularity). 

Proof : left as an exercise. □ 

Notice that the "circularity" condition is essential to have property ii) i.e. Tr[AiA2A3] = 
Tr[A2A3Ai], but Tt[AiA2A3] 7^ Tt[A2AiA3] 

Partial traces Rg e L{Hb) E L{Ha) of an operator R in L{Hi (E) H2) are 
defined accordingly as 

Rb = Tr^ [R] = ^^(<^„|i?|<^„),i Ra = Trs [R] = B(y;„|i? \ipn) b 



and circularity holds only for single-system operators, e.g., if R2 G L{Ha ® Hb), 
AeL(HA), B € L{Hb) 

Tr^ [A®lRiR2\ = ^a„(a„|i?ii?2|a„) = Tr^ [RiR2A®l] 

n 

TiA [A(EBRiR2\ = ^a„(a„|I(8)Bi?ii?2|a„) = Tr^ [1® BRiR2A®l\ 

n 

7^ a„(a„|i?ii?2 I ® S|a„) = Tr^ [R1R2 A®B] 



Exercise 13 Consider a generic mixed state g G L{H (g) H) and write the matrix 
elements of the two partial traces in terms of the matrix elements of g. 

Exercise 14 Prove that also partial trace is invariant under change of basis. 



B Uncertainty relations 

Two non commuting observables [X, y] 7^ do not admit common eigenvectors and 
thus it not possible to find common eigenprojectors and define a joint observable. 
Two non commuting observables are said to be incompatible or complementary, since 
they cannot assume definite values simultaneously. A striking consequence of this 
fact is that when we measure an observable X the precision of the measurement, as 
quantified by the variance {AX'^) = {X'^) — (X)^ , is influenced by the variance of any 
observable which is non commuting with X and cannot be made arbitrarily small. 
In order to determine the relationship between the variances of two noncommuting 
observables, one of which is measured on a given state {ip), let us consider the two 
vectors 

= {X - {X)m \i:2) = {Y - {Y)M , 
and write explicitly the Schwartz inequality (V'i|V'i)(V-'2|'02) > |(V'i|V'2)|^j i-e. [H] 



1 r 



|2 , i/^\|2 



{AX'){AY')>- \{F)\^ + \{C)r >-\{C)\\ (20) 



1 



|2 



4 



where [X, Y] = iC and F = XY-YX-2{X) (Y) . Ineq. represents the uncertainty 
relation for the non commuting observables X and Y and it is usually presented in 
the form involving the second inequality. Uncertainty relations set a lower bound to 
the measured variance in the measurement of a single observable, say X, on a state 
with a fixed, intrinsic, variance of the complementary observable Y (see Section 13.31 
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for the relationship between the variance of two non commuting observables in a 
joint measurement). The uncertainty product is minimum when the two vectors {ipi) 
and |'02) are parallel in the Hilbert space, i.e. l-^i) = —iX\'ip2) where A is a complex 
number. Minimum uncertainty states (MUS) for the pair of observables X, Y are thus 
the states satisfying 

{X + i\Y)\i:) = {{X)+tX{Y))\i,). 

If A is real then (F) = 0, i.e. the quantities X and Y are uncorrelated when the 
physical system is prepared in the state If |A| = 1 then {AX-^) = {AY"^) and 
the corresponding states are referred to as equal variance MUS. Coherent states of a 
single-mode radiation field |15j are equal variance MUS, e. g. for the pair of quadrature 
operators defined hy Q = -h={a^ + a) and P — -j=[a^ — a). 



